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HEAT TRANSFER IN THE TURBULENT INCOMPRESSIBLE BOUNDARY LAYER 

II - STEP WALL- TEMPERATURE DISTRIBUTION 
in 

2 . By W. C. Reynolds, W. M. Kays, and S. <7. Kline 

SUMMARY 

Heat- transfer rates and temperature profiles for the turbulent in- 
compressible flow of air over a flat plate with a stepwise temperature 
distribution (unheated starting length) were measured for a variety of 
step positions at Reynolds numbers up to 3.5 xiO^, Comparison of the data 
with existing heat- transfer analyses indicates that an improved analysis 
is needed. An integral analysis is made that agrees very well with the 
data and allows a simple correction for the unheated starting length. 

In addition, a differential analysis is made that allows prediction of 
the temperature profiles from the velocity profiles, and good agreement 
with experimental profiles is obtained. 


INTRODUCTION 

This report is the second of a series of four covering a three-year 
investigation of heat transfer in the turbulent incompressible boundary 
layer with arbitrary wall temperature (see ref. l) . The first report 
describes the experimental apparatus and presents results of experiments 
with constant wall temperature (ref. 2). Results of experiments and 
analyses for a step temperature distribution are presented herein. In 
part III the step-function analysis is used to predict heat- transfer 
rates for several cases of variable wall temperature, and the predic- 
tions are compared with experimental data (ref. 3). A simple method for 
handling arbitrary wall- temperature problems is presented. Part IV pre- 
sents an analysis of the effect of the location of the transition point 
on heat transfer in the turbulent boundary layer and compares the re- 
sults with experimental data (ref. 4) . 

The problem of heat transfer from a flat plate with a stepwise tem- 
perature distribution is of interest for two reasons. The "unheated 
starting length" problem is in itself of interest; but more important, 
the step temperature distribution provides the basis for analysis of 
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problems involving more complex arbitrary wall- temperature distributions 
by use of superposition (ref. 3) . Accurate analyses of the problem of 
step temperature distribution and experimental confirmation of these 
analyses are therefore of considerable importance. Prior to the present 
investigation, the only measurements of local heat-transfer coefficients 
with a step temperature distribution were those of Scesa (ref. 5). These 
data compare favorably with the analyses of Seban (ref. 5) and Rubesin 
(ref. 6) over the range attainable with Scesa 's experimental apparatus. 
However, this range was limited by the size of his test plate. The data 
presented in the present report agree with the existing analyses in the 
range of Scesa 's data but depart significantly from the analyses else- 
where. The largest departure occurs close to the discontinuity in the 
wall temperature, and this is the most critical region, with respect to 
both departure from isothermal performance and use in superposition solu- 
tions of more general problems. Thus an improved analysis is very 
desirable. 

The integral analysis presented in this report represents an im- 
provement over the former analyses in that the number of approximations 
has been reduced. The analysis involves assumption of the forms of the 
velocity and temperature profiles and an assumption about the mechanism 
of heat transfer. Use of the energy integral equation then leads to a 
differential equation for the thermal boundary- layer thic kn ess which, 
upon solution, allows calculation of the heat-transfer rates. This 
analysis is in excellent agreement with both the present data and all 
other previous data known. 

As just noted, it is necessary in an integral analysis to assume 
the shape of both the velocity and temperature profiles. In the present 
analysis, as in most such analyses, the two profiles have been assumed 
to be similar in shape. However, the velocity and temperature profiles 
in the boundary layer just downstream of a step in the wall temperature 
were carefully measured in the present investigation and were found to 
be dissimilar in this region. If the velocity profile alone is assumed, 
the differential energy equation of the boundary layer may then be 
be solved in detail. This solution yields the temperature profile of 
the boundary layer downstream of the step in wall temperature. This 
analysis is referred to as the differential aralysis and is presented in 
detail in this report. The results of the differential analysis and the 
measured temperature profiles agree well, and the heat-transfer predic- 
tions are in reasonable agreement with the experimental data. 

The present integral and differential anelyses are restricted to 
fluids having Prandtl numbers of unity, but should be adequate for 
fluids with Prandtl numbers near unity. Recently Ferrari made a simi- 
lar differential analysis in which he included the effect of Prandtl 
number and compressibility (ref. 7). This represents a considerable 
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extension beyond the present analyses but is extremely complicated. 
However, if the flow is incompressible, Ferrari ! s result reduces to the 
result of the integral analysis for any Prandtl number. It therefore 
appears that the assumption of Prandtl number of unity does not seriously 
limit the analyses presented herein. It is felt that the present analy- 
ses are useful because of their simplicity and that they adequately repre- 
sent the heat-transfer rates and temperature profiles. 

An essential feature of the integral analysis is a unique method of 
handling the shear stress and heat flux at the wall. Although profiles 
having infinite gradients at the wall are employed, turbulent viscosities 
and conductivities are constructed so that the wall shear and heat flux 
are finite and have their correct values. The idea of using a velocity 
profile with an infinite gradient at the wall is also employed in the 
differential analysis. This technique eliminates the need for separate 
consideration of a laminar sublayer and greatly simplifies the analyses. 

This investigation was carried out at Stanford University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 
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p fluid density, lb/cu ft 

a variable of integration 

T bJ shear stress in boundary layer, lb/sq ft 

shear stress at wall, lb/sq ft 


AMLYSIS 

Summary of Analyses 


The problem of turbulent heat transfer from a flat plate with a 
step temperature distribution (unheated starting length) has been treated 


approximately in a number of ways, 
general, be put in the form 

St 

Stiji 


The results of the analyses can, in 



where Stp represents the local Stanton number for heat transfer from a 
plate at constant temperature* The function f(z/x) thus represents a 
correction that may be applied to any suitable expression for the heat 
transfer from an isothermal surface* Reference 2 shows that the heat- 
transfer - momentum- transfer analogy of von Karmen (ref* 8) may be com- 
bined with the friction analysis of Schultz -Grunow (ref* 9) to give 


St 


T = 


1.60(ln Re x ) -2 * 58 

1 + 1.26(tn Re x )" 1,29 [5Pr + 5 ln(5Pr + l) - 14] 


( 2 ) 


This expression is felt to be the best available at the present t im e for 
turbulent incompressible flat-plate heat transfer. For air, having 
Prandtl number 0.7, reference 2 shows that equation (2) may be repre- 
sented in the range 10 8 < Re x < 10^ by 


St^Pr 0 * 4 = 0.0296 Re~° 


•fj 


■ 0.4 


( 3 ) 


Equation (3) represents a modification of the familiar Colburn analogy 
(ref. 10) . In using equation (3), the fluid properties appearing in the 
Stanton, Prandtl, and Reynolds numbers are to be evaluated at the free- 

stream static temperature; the factor (T w /Tj-°* 4 


corrects for the 





4995 


7 


effects of temperature-dependent fluid properties. Equation (3) is much 
simpler than the von Karman result (2) , and therefore is useful for com- 
putation purposes, 

A first approximation to the step temperature distribution problem 
can be obtained if the thermal effects are assumed independent of the 
hydraulic effects, in which case the heat-transfer rate can be determined 
directly from equation (3) by forming the Reynolds number with the dis- 
tance x - l instead of with x. This results in 


ft-Msr “> 

A second approximation was proposed by Rubesin (ref. 6), who at- 
tempted to separate the thermal and hydraulic effects in a semiempirical 
manner. Rubesin assumed l/7-power velocity and temperature profiles 
and, using the energy integral equation of the boundary layer, arrived 
at a differential equation between 6rp and x which, strictly speaking, 

applies only if the plate is at constant temperature. His differential 
equation contained a parameter m with which he proposed to separate 
thermal and hydraulic effects. In comparing his result with the limited 
data of Scesa (ref. 5), he chose the parameter m to be zero. His 
final result may be put in the form 


% 1 - erf* 

Seban (ref. 5) used a l/7-power velocity profile and a temperature 
profile that was linear near the wall. He "patched" the linear portion 
to the l/7-power profile in the outer portion of the boundary layer. 
Then, by using the energy integral equation, Seban arrived at the fol- 
lowing result: 


= p r - 2 / 9 [i . m 9 / lo r 1/9 

st T L w J 


( 6 ) 


Except for the Prandtl number dependency, this result is identical with 
the integral analysis presented in the present work. One would expect 
that, well downstream of the step, the effect of the step should die out, 
and St/Stij should approach unity. Thus the factor Pr - 2/9 gives an 
improper limiting behavior for x -► •», This erroneous Prandtl number 
dependency in Seban* s solution is a result of his "patching" of the tem- 
perature profile. 
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In using integral methods, there is always a problem as to what to 
use for the value of the heat flux at the wall (see eq. (l6)). Rubesin 
in effect modified the isothermal value, while Seban used the slope of 
the temperature profile in his laminar sublayer. The latter approach is 
certainly more desirable, but leads to the introduction of an incorrect 
dependency on Prandtl number. In the present integral analysis, this 
term is evaluated by a limiting procedure as the product of the slope of 
a power profile (which is infinite at the wall) with an eddy diffusivity 
(which is made zero at the wall in such a manner that the shear stress 
and heat flux have their correct values) . Then, after assuming l/m- 
power velocity and temperature profiles, one obtains 


St 

Stiji 



(m+2) /5(m+l) 


f 


-l/ (m+2) 


(?) 


For m = 7 this result is similar to Seban* s but does not have the 
erroneous Prandtl number dependency. Equation (7) with m = 7 is in 
excellent agreement with the present data, and is recommended as the 
best step temperature distribution analysis. 

The chief objection to the integral methods is that the forms of 
the velocity and temperature profiles must be assumed to be similar. It 
is well known that this is approximately true for a constant- temperature 
plate, but it is not true for nonisothermal heat transfer. To avoid this 
assumption it is necessary to use the energy equation in differential form 
and to solve directly for the temperature profiles. This has been done in 
the present differential analysis. The heat-transfer rates may then be 
determined from the slope at the wall, and the result may be written as 



where C n and X n are expansion coefficients and eigenvalues, re- 
spectively. The differential analysis is not in as good agreement with 
the experimental data as the integral analysis, but it does allow pre- 
diction of the temperature profiles and thus is of interest. 

Recently Ferrari (ref. 7) has made a similar differential analysis, 
except that he allowed for a sublayer. Ferrari included compressibility 
effects, fluid-property variations, and Prandtl number effects, and thus 
his is probably the most advanced analysis available at the present time. 
However, for incompressible flow, Ferrari * s result reduces to the result 
of the present integral analysis (eq. (7)), which is obtained far more 
easily. 
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In addition to the several analyses, two empirical correction formu- 
las have been proposed. Maisel and Sherwood (ref. ll) made mass-transfer 
measurements from a flat plate, and Klein and Tribus (ref. 12) examined 
their data and proposed the empirical relation 



( 9 ) 


Note that this equation is quite similar to the result of the present 
integral analysis. Jacob and Dow (ref. 13) made mean heat-transf er meas- 
urements on a cylinder in axial flow and found that the correlation of 
their mean heat-transfer data implied the local Stanton number could be 
given by 


St 

st T 


/7 

\ / *T 

\2.75 /i 

r\3.75 

0.8 + 0.2(± 

00 

r- 

• 

O 

i 

) + ia8 U 

0 


( 10 ) 


These empirical equations are in fair agreement with the present data. 


Integral Analysis 

The momentum equation for the incompressible boundary layer on a 
flat plate may be written in the following form (see appendix D of ref. l) 




Pu £ 


t (i - 1) ■ * + ^r © 2 * - tH 3 (t) * 


( 11 ) 


Note that for y » 0 equation (ll) reduces to the more familiar momentum 
integral equation 


2 


w 


pu; 




(lla) 


Moreover, if y is set equal to 5 in equation (ll), the correct result, 
T-^-j(x,&) = 0, is obtained. Now if it is assumed that the velocity profile 

may be represented by an expression of the form 


- = te 




m 


u 

u 


( 12 ) 
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the shear-stress distribution in the boundary Layer may be determined 
from equations (ll) and (lla) as 


/y\(2-Hn)/m 

“ 1 " UJ 


(13) 


The mechanism for momentum diffusion in the turbulent boundary layer may 
be described by 



P 



(14) 


Strictly speaking,, equation (l4) holds only far away from the wall where 
the viscous effects are negligible and all the momentum transfer is by 
turbulent eddies. However, the power profile assumed by equation (l2) 
will have an infinite gradient at the wall, since m >1; if laminar 
effects were included near the wall, the shear stress would have to be 
infinite. It is therefore necessary to neglect the laminar terms and 
construct an eddy diffusivity such that the wall shear stress has its 
correct value. This is done simply by combining equations (l2) , (l3), 
and (14) : 


Lm C f & 

v - m — - Re x 



(15) 


Note that the eddy diffusivity is zero both at the wall and at the outer 
edge of the hydraulic boundary layer. 

The energy equation for the boundary layer may be written as fol- 
lows if dissipation is neglected and the temperature difference t - t 

is assumed constant in the region of interest: 



(1 - 0) dy 


(16) 


The mechanism for turbulent heat diffusion in the boundary layer may be 
described by 



pc 


P 


st bi 

- e H “57* 


(17) 


It will now be postulated that the eddy diffusivity for heat is equal to 
the eddy diffusivity for momentum, which is the familiar Reynolds analogy. 
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But, since equations (14) and (l7) are assumed to hold everywhere in the 
boundary layer, even near the wall, this postulation implies that the 
Prandtl number is 1. This is no serious restriction on the analysis, 
since the result will simply be a correction on the isothermal Stanton 
number, and the Prandtl number dependence can be later incorporated in 
the expression for the isothermal Stanton number (3) • The final assump- 
tion of this analysis is that the temperature profile is similar to the 
velocity profile, based on its own boundary-layer thickness: 

( ' \l/m 

fi) < ie > 

Note that no other power profile could be used, or the heat flux at the 
wall would be infinite. Combining (15) and (l8) with (17) gives 


PC p U » U w " U 


1 - 


(24m) /] 


m 


-1/. 


m 


(19) 


Setting y = 0 in equation (19) gives 

..4 ft)'* 

This expression will now be used as the left side of equation (l6) • 
Then, for the assumed velocity and temperature profiles. 


( 20 ) 


foy 1/m - (- 

& T \( 1 / m )+i 

V 5 / V 

~) 


Integration of (2l) yields 

^(x) ^> T /S 

dS 
5 

/S(Z) OO 


i - twb If 


1 /“ d(o. r /b) 


dx 


( 21 ) 


m 


1 


-m 


2/m 

6rp\(2/ m )+l 


m+1 ^ a Q?) 


( 22 ) 


where S(Z) is the hydraulic boundary- layer thickness at x = l, 
integrations result in 


The 


6 ( x ) f /iy.( 2 *)/ml-( m+1 >/( m+2 > 

6TTT ' L 1 ' Uv 


T 


(23) 
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which is equivalent to 


8kp 

T 




(m+2) / (m+l)) ■*l/(2+m) 

) 


(23a) 


It is well known that in the range 10^ < Re x < 10^ the boundary- layer 

4/5 A 

thickness 6 varies as x ' , and thus equation (23a) may be written as 


5 


-(I) 


4 (m+2) /5(m+l) 


-l/(2-nn) 


(24) 


Substituting (24) in (20) gives 

st - T |_! - (i) J < 25 > 

Note that, far away from the discontinuity in wall temperature, where 
l/x -*• 0, St -> Strp = Cj/2. Thus, the final result may be written as 


St 

St, 


T 




4(m+2)/5(m+l) 


•l/(2+m) 


(26) 


If m is taken as 7, equation (26) reduces to 

1-1/9 


St 

Stn 


er 


(26a) 


The mean Stanton number at the end of a plate of length L is de- 
fined as the mean St over the heated portion of the plate } 


St m l 


r 

- I Jr 


St(x) dx 


(27a) 


From equations (26) and (3) (for m = 7), 


•L Jl 


^ - zh ( I )’ 0 ' 2 / °' 0 - 2 (* - ‘- 9/l0 )" l/9 - (™> 


S66^ 
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The integral is readily evaluated by setting o 


9/l0 


z , and finally 


Stiji 



(28) 


Note that, as 


l/i 


Stm ^ 5 
StT 1 


(28a) 


which is the well-known result for the constant- temperature plate (ref. 
14=). In equations (27b) to (28a) , St T is to be evaluated at L. 


Differential Analysis 

The energy equation for the turbulent incompressible boundary layer 
may be written, neglecting dissipation and assuming that t w - is 

constant in the range of interest, as 


30 30 3 ( OQ \ 

u s + v SJ * C29) 

Again it is assumed that the heat- transfer mechanism everywhere in the 
"boundary layer is given "by equation (17). In this analysis it will again 
be assumed that the velocity profile may be given by an expression of 
the form 

(m * 5 to 8) (12) 

00 x 

It should be observed that, although the velocity gradient is infinite 
at the wall, a reasonable answer will be obtained. The continuity equa- 
tion for incompressible flow is 


du dv 
dx * ciy 


0 


Combination of equations (12) and (30) gives 

_L = _1_ ^ /y\( 1 / m ) +1 
u_ ~ m+1 dx \5/ 

OO 


(30) 


( 31 ) 
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The eddy diffusivity for heat was calculated for the integral analysis 
and is^ assuming e H = e M-> 


C H C f 5 f /y\( 24ni )/ m 

V = m — X Re x 1 - 

Note that the assumption = everywhere in the boundary layer , in' 

eluding the region near the wall, again implies Pr = 1. Now, if new 
coordinates 



K = x/z 


n = y/8 

are introduced, equation (29) transforms to 


= /x\ i *h dfe /x\ i be Se H 1 

u . ^ u^ 8 dx 1 ^ u_ 6 drj \8 j 8 ^ 2 + \6j 8 dr^ u 

(32) 

Substitution of (12) in (lla) results in 


= m d8 d8 

2 (m+1 ) (m+2 ) dx dx 


(33) 


where M = m/(m+l) (m+2) . Blasius has proposed that, for turbulent flow 
over a flat plate, the friction factor and boundary- layer thickness may 
be related by (see ref. 14) 


Cf 

— = 0.0228 Re 



(34) 


which, when combined with (33), leads to 

| = (0.0228) 4 / 5 (0.8)“ 4 / 5 k" 4 / 5 Re"°* 2 


(35) 


and 

= (0.0228) 4 / 5 (0. 8) ^K^Re" 0 * 2 


(36) 
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Substituting (l2), (3l) , (15), (35), and (36) into (32) gives 

rji/m^ + o.8M(m-l)^ 1 / ra ) +1 - b/m) - TJ 2+( 1 /m)J |£| 

(37) 


With the change of coordinates 


w = Zn £ 


= n = 

u 


equation (37) can be reduced to 

S0 _ 0.8M 


.2 

/_-m r,2) ^ 0 

^ ~ U j ^2 


(38) 


This linear parabolic differential equation resembles the equation de- 
scribing transient heat conduction in a slab, and techniques similar to 
those used in such problems will be used in its solution. Equation (38) 
is to be solved in the region w > 0 for 0 < z <1. The boundary 
conditions for 0(w,z) are 


6>(w,0) = 0 
e(w,i) = i 


(39) 


In addition, there is the "initial” condition that the boundary- layer 
air must be at the free-stream temperature at the start of the heated 
section; that is, 


0(0, z) = 1 


(40) 


In handling transient conduction problems it is convenient to divide 
the solution into "transient” and ”steady-state" components* The same 
technique is useful here. The "steady-state” solution is simply the tem- 
perature profile far downstream of the temperature discontinuity, 



(41) 


V 


Note that 0 = z satisfies equation (38) 
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Thus, 6 is divided into 


0 = z + <P (w,z) 


(42) 


where <p(w,z) is to satisfy 


d<P 0_ 1 8M / m _ 2 ) 

m dz 2 

subject to the boundary conditions 

«P(w,0) = 0 

«p(w,l) = 0 

and the "initial" condition 

<p(0,z) = 1 - z 

Equation (43) is solved by separation of variables, where it is 
assumed that 


(43) 


(44) 


(45) 


<p = W(w)z(z) 


which leads to 


Sf - ^ (z - m _ *2) !£ . „ X 

The resulting differential equation for W(w) is 

W = -XW 

which has the solution 

W = Ce 

The differential equation for Z(z) is 


-Xw 


Z" + 


m 


0.8M( 


z _m - z 2 ) 


XZ = 0 


(46) 


(47) 


(48) 


(49) 


(50) 


The transformation £ = z m+2 reduces equation (50) to 

d 2 Z /m+l\ f . M dZ (m+1) X 

d ^2 + (m+2/ ^ ^ d£ + (5+2) 0.8 Z “ 


( 51 ) 
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which has the solutions (ref. 15, p. 178) 


where 


to 

i 

r* 

o 


z = c 1 F(a,p;r;U + 

c/'b 

a “ 2 (m+2) L 1 

- vTT 

B - - 1,1 

1 + y/T 

2 (m+2) 1 


r = 


m+l 

m+2 


(52) 


Now from (44) , 


Z(0) = 0 


and thus 


C 


1 


= 0 


From (44), 


z(i) = 0 = c 2 F(a-r+i,0-r+ij2-rji) 

But it is well known that (ref. 16) 

\ r (c)T (c-a-b) 

F(a/b,c;l) - r ( c _ a )r(c-b) 

and thus, from (53), since Cg is not zero. 


(53) 


(54) 


0 = 


r(£i±h rm 

\m+2/ 

/ 2m+5-B \ 7 2m+5+B \ 

\ 2m+2 / \ 2m+2 ) 


(55) 


where 


B = y/l + 5X (m+l) (m+2) 
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Since the gamma function of a negative integer is infinite, the 
eigenvalues are given by 


2m + 5 - B 
2m + 2 


= -n + 1 


n = 1,2,. 


and thus the eigenvalues are 


X - [2(n-l)(m+2) + 2m + 5 ] 2 
5(m+l) (m+2) 


- 1 


n 


The solution for <p may now be reduced to 


^n e ^ nW zF(n-l+a,-n;a;z m+2 ) 

n=l 


( 56 ) 


(57) 


(58) 


where a = (m+3)/(m+2). 

It is now necessary to expand 9 as a series of eigenfunctions in 
order to satisfy the initial condition (45) . It may be shown that 
( re T. 15, p. 225), if Zj and are solutions of (50) associated 

with the eigenvalues and \, respectively, the solutions Zj and 

have the following orthogonality property: 

2 z j z k dz = o :-f j i k (59) 

“ z 

Thus by multiplying (58) by 

g zF(n-l4a,-n;a;z n + 2 ) 

Z ” z 



setting w = 0, and integrating from z = 0 to z = 1, the constants 
C n may be evaluated as 


C 


n 



zF(n-l+a,-n;a;z m+2 ) dz 



z 2 F 2 (n-l+a,-n;a;z g ' +2 ) 
z _gl - z 2 


(60) 


dz 
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CD 

CD 


Since the second argument of the hypergeometric function is a negative 
integer, the functions are actually polynomials, and the coefficients may 
be evaluated by expansion and integration term by term. This has been 
done for n = 1, 2, and 3, resulting in 

C-^ = (a - l) (a + l) 


C 


2 


a 2 (a - l) (a + 3) 
4 


C 


3 


a^(a - l) (a + l)^(a + 5) 
36 


Evaluation of the higher order coefficients is tedious and is best 
accomplished by numerical methods. This has been done for m = 5.6, 
which applies in the Reynolds number range of the present data, and it 
was found that 


C 4 = 0.1110 
C 5 = 0.0940 


The final solution for 0(x/Z,y/S) may be written as 


where 


e = ( 

\l/m 

i) 1 

ft / 

- m 

i + 2 cM 


;)' 


° J 

n=l 



G n(i) » F n-l-fa,-n;aj(|) (m+2)/m 


(61) 


(62) 


The heat- transfer rate may be found by differentiation of (6l) and 
use of (17) at y = 0. This results in 



(63) 


The hypergeometric functions of interest in determining the first 
five eigenfunctions for m » 5.6 have been calculated to four decimal 
places with an IBM 650 digital computer. These functions are given in 
table I and are shown in figure 1. 
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RESULTS AND DISCUSSION 

Comparison of Analyses with Heat-Transfer Experiments 

A heated plate with an active flow length of ahout 5 feet was used 
in this investigation. The heated surface cor. sisted of 24 individually 
heated strips, thermally insulated from each other. Within control 
limitations, an,y desired surface temperature or heat-flux distribution 
could be set up and desired heat- transfer measurements obtained. The 
experimental apparatus and techniques are described in detail in refer- 


Eighteen runs were made for which a step wall temperature existed. 
Data were obtained at flow Reynolds numbers up to 3.5 X 10^, with step 
Reynolds numbers ranging from 0.5 X 10^ to 3 X 10^ . These data are pre- 
sented in table Il(a). The discontinuities were quite sharp, and ex- 
cellent "steps" were obtained. In reducing these data, allowance was 
made for conduction between the strips within the plate, and this cor- 
rection was only important for the first heated strip. The experimental 
Stanton numbers are actually the average Stanton number over each strip, 
since the strips are of finite flow length. An analysis indicated that 
the Stanton number at the center of the strip, where the Reynolds number 
is calculated, does not differ significantly from the average over the 
strip, the largest difference occurring on the first heated strip, where 
the average Stanton number exceeds the Stanton number at the strip center 
by only about 2 percent. For all other strips, the difference is less 
than l/2 percent. Analysis indicated that the imperfectness of the dis- 
continuity in the surface temperature reduced the heat transfer on the 
first heated strip from the abrupt-discontinuity value by about 2 per- 
cent and had no appreciable effect on the heat transfer from the strips 
farther downstream. Since the error in averaging was approximately 
compensated for by the error due to the step imperfection, no correc- 
tion was made for either. 

The data are shown graphically in figure 2. It can be shown (ref. 
17) that an adequate method for considering the influence of temperature- 
dependent fluid properties for gas flow in both external and internal 
boundary layers is to evaluate all fluid properties at the free-stream 
static temperature and then to include all temperature-dependent property 
effects in a factor (T w /T^) m , where the exponent m is a function of 

geometry alone. Examination of the recent results of an analysis of 
Deissler and Loeffler (ref. 18) indicates that, for the turbulent in- 
compressible boundary layer, the Stanton number varies as (T w /T )"^'^, 

OD 

other things being equal. Thus the results are presented in the form 
St(T w /T g* 4 against Re x . 
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21 


Reference 2 shows that heat transfer to air from isothermal plates 
may be represented by 


StrpPr 


0.4 


= 0.0296 Re” 0 ’ 2 



( 3 ) 


Equation (3) is plotted in figure 2. The effect of the unheated start- 
ing length is considerable (especially near the step, where the data 
show a correction as large as 50 percent of St j) , and the data gradually 

approach the isothermal correlation downstream of the wall- temperature 
discontinuity. 


The data are compared with the various analyses in figure 3. The 
ratio of the local Stanton number to the isothermal Stanton number as 
given by equation (3) is plotted against the factor 1 - (l/x) a / 10 , 
which alloys equation (26a) to be plotted as a straight line with a -l/9 
slope. The simple analysis that neglects the dependence of heat transfer 
on the hydraulic boundary layer (4) is high near the step, as is the 
semiempirical analysis of Rubesin (eq. (5)). The Seban result (eq. (6)) 
lies parallel to the data but is displaced by the Pr “2/9 factor. The 
integral analysis employing l/7-power profiles (eq. (26a)) is in excel- 
lent agreement with the data over the entire range. However, in the 
Reynolds number range of these tests, the velocity profile is actually 
closer to a 1/5. 6-power profile (see fig. 7). The differential analysis 
for m = 5.6 is slightly high near the step because, near the step, 
most of the thermal effects are confined to the region near the wall, 
where the assumed velocity profile is least accurate (see fig. 9(a)). 
Downstream, where the power profile is a good approximation, the differ- 
ential analysis for m = 5.6 is in good agreement with the data. The 
success of the integral analysis for m = 7 probably lies in the fact 
that the temperature profiles just downstream of * the step are consider- 
ably "fatter" than the velocity profiles, and thus are actually closer 
to l/7-power profiles (see figs. 8 and 9). 

Since the agreement of the present integral result (eq. (26a)) with 
the experimental data is so good, equation (26a) is recommended as the 
best correlation available at the present time for a step temperature 
distribution in turbulent incompressible flow. The data are shown as 
correlated by this equation in figure 4. The correlation is excellent 
over the entire test range. 

Hie data of Scesa (ref. 5) have been correlated in a similar manner, 
and this correlation is shown in figure 5. Most of Scesa T s data were 
obtained for a double step in the wall temperature, and thus a more com- 
plicated correction was necessary. For this reason the "corrected 
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Nusselt number" is plotted against Re x , the correction method following 

the double-step example of reference 3. These data were obtained over a 
rather limited range, as is indicated by figure 5, and there is con- 
siderable scatter in the data. The accuracy of data obtained from the 
present apparatus is believed to be considerably better than that ob- 
tained with Scesa's small plate. This view is substantiated by the 
presence of less scatter in the present date. 

The mean heat- transfer measurements of Jacob and Dow (ref. 13) are 
compared with the present analysis (eq. (28)) in figure 6. These data 
were obtained from tests on a cylinder in axial flow, and the Reynolds 
number was low so that transition effects may have been of some importance. 
In consideration of these facts, the agreement of these data with equa- 
tion (28) is believed to be very satisfactory. These data are also pre- 
sented on the basis of a corrected mean Nusselt number against Rex. 


Comparison of Survey Data with Theoretical Results 

Figures 7 and 8 show velocity and temperature profiles obtained for 
constant wall temperature, as reported in reference 2. The profiles are 
similar and may be represented very well by a l/5. 6-power formula. The 
apparatus used to obtain these data and the survey data presented here- 
in is described in detail in reference 2. 

Temperature surveys with a step wall-temperature distribution were 
obtained at three points on the plate and at two different velocities 
for each point. In these tests the first 11 strips were unheated, and 
the remaining strips were held at constant temperature. The data from 
these surveys are shown in table II (b) and are compared with the pro- 
files predicted by the differential analysis (m = 5.6) in figure 9. In 
these comparisons the thermal boundary- layer thicknesses were determined 
in the manner described in reference 2; the conduction thickness 6* 

was determined by numerical integration of the temperature profile, and 
the "best power fit" to the data was determined by plotting 9 against 
y on log- log paper. The conduction thickness is related to the 

thermal boundary- layer thickness for a power profile of the form 
e = (y/Brp) 1 / m by the equation 

5rp = (l + m) S^i 

The values of Sqi determined in this manner were quite close to the 
99-percent thicknesses (ref. 2). In determining the thermal boundary- 
layer thickness from the differential solution, the same techniques 
were applied. 





23 


The surveys closest to the step (x/l = 1.20) show why the heat- 
transfer rates predicted from the differential analysis are not too good 
near the step. Wear the step, most of the thermal effects are confined 
to the innermost portions of the boundary layer, where the velocity 
profile assumed for the analysis is least accurate. The predicted and 
experimental profiles near the step are not in agreement, as is shown 
by figure 9(a) . Farther downstream the predictions are remarkably ac- 
curate, as shown by figures 9(b) and (c). These survey data also show 
why the integral analysis employing 1/7-power profiles works so well; 
the temperature profiles downstream of the temperature discontinuity are 
much fatter than the velocity profile, and in fact are quite close to 
being l/7-power profiles. This is probably responsible for the success 
of the integral analysis with m = 7. 

CONCLUDING REMAEKS 

An integral analysis has been obtained that allows calculation of 
heat transfer from a flat plate with a step wall -temperature distribution 
(eq. (26a)). The analysis is in excellent agreement with experimental 
heat-transf er data of the present and earlier investigations. A differ- 
ential analysis was obtained that allows the prediction of the tempera- 
ture profiles in the boundary layer downstream of the temperature dis- 
continuity. The predicted profiles are in good agreement with experi- 
mental profiles, but the heat-transf er rates predicted by this analysis 
are slightly high near the temperature step. 

The method of handling the heat-transf er rate at the wall used in 
the integral analysis is unique in that it allows a power profile having 
an infinite gradient at the wall to be used everywhere in the boundary 
layer, while still yielding good results. This eliminates cumbersome 
’patching" of profiles to obtain a finite slope at the wall. The same 
technique was used in the differential analysis, wherein the energy 
equation was solved in spite of the infinite velocity and temperature 
gradients at the wall. It is believed that this method may also be use- 
ful in the treatment of more complex problems, such as flows with pres- 
sure gradients, since it provides both analytical simplicity and accurate 
prediction in the cases thus far explored. 

The recommended correlation for heat transfer in the case of a step 
wall temperature is the result of the integral analysis, 

& 4 - erf' 


(26a) 
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The isothermal Stanton number St|p for incompressible flow may be com- 
puted from the relation 


Strp 


= 0.0296 Re" 0,2 Pr" 0 * 4 

X 



( 3 ) 


In this relation the fluid properties are to be evaluated at the free- 
stream static temperature. 

t; 


Stanford University, 

Stanford, Calif., October 22, 1957. 


REFERENCES 

1. Reynolds, W. C.: Heat Transfer in the Turbulent Incompressible 

Boundary Layer with Constant and Variable Wall Temperature. Ph.D. 
Thesis, Stanford Univ., 1957. 

2. Reynolds, W. C., Kays, W. M., and Kline, S. J.: Heat Transfer in the 

Turbulent Incompressible Boundary Layer. I - Constant Wall Temper- 
ature. NASA MEMO 12-1-58W, 1950. 

3. Reynolds, W. C., Kays, W. M., and Kline, S. J.: Heat Transfer in the 

Turbulent Incompressible Boundary Layer. Ill - Arbitrary Wall Tem- 
perature and Heat Flux. NASA MEMO 12-5-58W, 1958. 

4. Reynolds, W. C., Kays, W. M., and Kline, S. J.: Heat Transfer in the 

Turbulent Incompressible Boundary Layer. IV - Effect of Location 
of Transition and Prediction of Heat Transfer in a Known Transition 
Region. NASA MEMO 12-4-58W, 1958. 

5. Scesa, S.: Experimental Investigation of Convective Heat Transfer to 

Air from a Flat Plate with a Stepwise Discontinuous Surface Temper- 
ature. M.S. Thesis, Univ. of Calif*, L951. 

6. Rubesin, Morris W.: The Effect of an Aroitrary Surface -Temperature 

Variation Along a Flat Plate on the Convective Heat Transfer in an 
Incompressible Turbulent Boundary Layer. NACA TN 2345, 1951. 

7. Ferrari, Carlo: Determination of the He it Transfer Properties of a 

Turbulent Boundary Layer in the Case of Supersonic Flow When the 
Temperature Distribution Along the Constraining Wall is Arbitrarily 
Assigned. Rep. No. CAL/CM-807, Cornell Aero. Lab., Inc., Mar. 1954. 
(Contract NOrd- 14523.) 





4995 


25 


i 

o 


o. von Karman, Th.: The Analogy Between Fluid Friction and Heat Transfer. 

ASMS Trans. , vol. 61, no. 8, Nov. 1939, pp. 705-710. 

9. Schultz-Grunow, F.: New Frictional Resistance Law for Smooth Plates. 

MCA TM 986, 1941. 

10. Colburn, Allan P.: A Method of Correlating Forced Convection Heat 

Transfer Data and a Comparison with Fluid Friction. Trans. Am. 

Inst. Chem. Eng., vol. XXIX, 1933, pp. 174-206. 

11. Maisel, D. S., and Sherwood, T. K. : Evaporation of Liquids into 

Turbulent Gas Streams. Chem. Eng. Prog., vol. 46, no. 3, Mar. 

1950, pp. 131-138. 

12. Klein, John, and Tribus, Myron: Forced Convection fron Nonisothermal 

Surfaces. Eng. Res. Inst., Univ. Michigan, Aug. 1952. (Contract 
AF 18(600) -51.) 

13. Jacob, Max, and Dow, W. M. : Heat Transfer from a Cylindrical Sur- 

face to Air in Parallel Flow with and without Unheated Starting 
Sections. Trans. ASME, vol. 68, no. 2, Feb. 1946, pp. 123-134. 

14. Eckert, E. R. G.: Introduction to the Transfer of Heat and Mass. 

McGraw r -Hill Book Co., Inc., 1950, p. 75. 

15. Hildebrand, F. B.: Advanced Calculus for Engineers. Prentice-Hall, 

1949. 

16. Erdelyi, A., Magnus, W. , Oberhettinger, F., and Tricomi, F. G.: 

Higher Transcendental Functions. Vol. 1. McGraw-Hill Book Co., 
Inc., 1953, p. 61. 

17. Kays, W. M. : A Summary of Experiments and Analyses for Gas Flow 

Heat Transfer and Friction in Circular Tubes. TR-22, Dept. Mech. 
Eng., Stanford Univ., June 30, 1954. (Contract NGonr 251.) 

18. Deissler, R. G., and Loeffler, A. L. : Turbulent Flow and Heat 

Transfer on a Flat Plate at j High Mach Number with Variable Fluid 
Properties. Paper No. 55-A-133, ASME, 1955. 



26 


TABLE I. - HYPERGEOMETRIC FUNCTIONS FOR 
STEP-FUNCTION ANALYSIS 


[F n (z) = F(n-l+a,-n;a;z) j a = (m+3)/(m+2); 
m = 5.6 . ] 


z 

F x (z) 

F^(z) 

F 3 ( z ) 

F 4 (z) 

f 5 ( 2 ) 

0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

.05 

.9500 

.8250 

.6240 

.3932 

.1571 

.10 

.9000 

.6509 

.3219 

.0011 

-.2318 

.15 

.8500 

.4972 

.0871 

-.2214 

-.3286 

.20 

.8000 

.3572 

-.0876 

-.3146 

-.2591 

.25 

.7500 

.2311 

-.2072 

-.5146 

-.1170 

.30 

.7000 

.1188 

-.2797 

- .3529 

.0322 

.35 

.6500 

.0204 

-.3114 

- . J569 

.1474 

.40 

.6000 

-.0642 

-.3087 

-.0490 

.2079 

.45 

.5500 

- o 1349 

-.2784 

.0526 

.2095 

.50 

.5000 

-.1919 

-.2269 

.1339 

.1609 

.55 

.4500 

-.2349 

-.1608 

.1859 

.0791 

.60 

.4000 

-.2642 

-.0868 

.2035 

-.0137 

.65 

.3500 

-.2796 

-.0115 

. 1865 

-.0946 

.70 

.3000 

-.2812 

.0586 

.1389 

-.1435 

.75 

.2500 

- .2689 

.1169 

.0692 

-.1472 

.80 

.2000 

-.2428 

.1567 

-.0095 

- . 1033 

.85 

.1500 

-.2028 

.1716 

-.0780 

-.0242 

.90 

.1000 

-.1491 

.1548 

-.:,202 

.0594 

.95 

.0500 

-.0745 

.0134 

-.0138 

.0123 

1.00 

0 

0 

0 

0 

0 
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TABLE II. - EXPERIMENTAL DATA SUMMARY 
(a) Heat- transfer data 



At„, Q.;. 


lb 

°F Bt u 

Btu 

( hr ) ( 3q ft) 
*10-" 

: < h r ) ( mt ft) 

(hrXfjq n)(Op) 


8 l 1 /rr y0.4 


Re, = 0.MS-10 6 ; . 


(hr)(sq ft) 
* 10-3 


Btu 


{ hr ) { »q ft) 


(hr) (oq ft )(°P) I 


0. 0.’34 i fc/cu ft. 









n , 







3 

25.4 

0 

___ 






’ !524 

34.0 

2.2 






4 

25.4 


— 


— 

— 

.445 

33.9 

21 .a 

553 

24-4 

3.00 

J . 04 

5 

25.4 

19.8 

394 

19.40 

3 . 20 

3 . 26 

. 569 

33.9 

22.5 

456 

£ . p 

2 . 48 

2 . L>2 

f* 

25.4 

19.5 

321 

16.10 

? . U h 

2.69 

.664 

33.7 

22.4 

423 

18.90 

2.33 

"? , * *' 

7 

25.4 

19.9 

304 

16.27 

2.50 

2.54 

. HO 5 

33 . 7 

22.2 

398 

17.93 

2.24 

2.27 

6 

25.4 

19.9 

295 

14.40 

2.3 7 

?. 4 0 

.922 

33.5 

22.3 

386 

1 7 . 30 

2.1S 

2.18 

9 

25.4 

20 . 2 

? h a 

1 4 * /i \) 

2.35 

2 . 3 9 

1 .045 

33.6 

22.3 

387 

17.34 

2.10 

2.19 

10 

25.5 

20.3 

262 

12. 87 

2.11 

2.14 

1 . 1 66 

33.5 

22.9 

360 

15 . 73 

1.96 

1 .98 

11 

25.5 

20 . 7> 

2 SO 

12.29 

? . 0 1 

2.04 

) .290 

33 . 5 

23 .0 

346 

15.05 

1 .87 

1 .90 

12 

25.5 

20 .4 

25? 

12.33 

2.01 

2.04 

1.411 

33 . 0 

6 

/, 

15.75 

1 . 9 6 


13 

25.5 

if j . 4 

252 

12.37 

2 . 02 

. Ob 

L . j2fi 

33.6 

) j j '* ■ ■! !i 

340 

15.05 

1 . 8 7 

1 . 67 

14 

25.4 

ffO .5 

252 

12.29 

2.01 

2.04 

1.656 

33.6 

2? . d 

345 

15.12 

1.86 

1 . 91 

15 

25 . 4 

20.5 

245 

11.92 

1.96 

1 . 99 

l . 761 

33 . 6 

23.0 

34 1 

1 4 . 63 

1 .85 


I’! 

25.5 

20 . U 

24 2 

11.72 

1 .92 

1.96 

l . 890 

;«,3 . s 

23 . : 

335 

14.50 

1.60 

1 .63 

17 

25 . 5 

20 . 6 

24? 

11 .73 

1.9? 

1 . ;H) 

2.01 

33 . b 

22.9 

329 

14.42 

1 . 80 

1 . 62 

1 H 

25.4 

20 . 7 ' 

2.34 

11.27 

1 . 8b 

1 . H K 

0. 12 

33 . s 

22.8 

322 

14.12 

1 . 76 

; 1.7H 

' 19 ' 

25.4 

20 . H 

23 5 

11.30 

1 . 65 

1 „ flh 

2 . 2 5 

3.5 . u 

23.0 

321 

1 13.95 

! 1 .74 


i 20 

25.4 

20 . d 

236 

11 . 13 

1.6?. 

1 . Mh 

? . 36 

33.4 

23 .0 

32b 

14.12 

1 . 76 

1 . 70 

21 

25.4 

20 . a 

229 

10.94 

l.auj 

; # f-j^i 

;■! . 4 8 

* 3.4 

123.4 

315 

13.47 

1 .66 


22 

25.4 : 

?: .o , 

210 

10.00 

'1.65 

i ► i..7 

2 . 6.0 

33.3 

23.5 

293 

12.47 i 

1 . 56 

1 *. 5-6 

23 

25.5 

21.2: 

220 

10 .39 

,l.h0| 

1 • 12 , 

2 . 7,3 

33.3 

23.4 

295 i 

12.60 j 

1 .58 j 

j 1 60 , 


R,1 l " 

0 , 509*10® ; l m = 

.•3. i b; Vm » u. 

. 0751 

ib/cu ft 


Re « 

l 

1.00*10®; t. = 84.2° F; p = O.C 

)736 : 

b/dli ft 

.. 

17.4 






1 ) 1 6 

or 






3 

17.4 

0 



mil 


1 

! 30,0 

- 





4 

17.4 

0 






30 .0 


0 




s 

1 7 . 4 


;-j 










6 : 

17.4 
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7 9*9 


X 




1 l 1 

17.4 

21 . L 

284 

13.4 7 

.5 . 23 

3 . 28 

.550 ; 

29. 9 

1 .9 

--- 







0. 259 
.429 
.555 
. 7 45 
• H95 

1 . 055 
1 . 200 
1 .560 
1.515 
1.670 

1 . 6 26 
1 .96? 
2.14 
2.25 
2.45 
2, PI 


5.07 

5.22 

3.57 


9 

17.4 

21.3 

246 

11.5? 

,2.76 

2. HO 

.714 

29.9 

20.8 

3 98 

19.14 

2 ' p 7 

:> ' 7 i 

1.207 

10 

17.4 

21.2 


10.19 

? . 44 

2.48 

. ?97 

29 . 8 

20.9 

34 8 

16.66 

2.33 

2.31 


11 

17.4 

21.3 

211 

9 . 8 7 

2 , 3 7 

2.41 

.678 

29 . 7 

PI.? 

331 

15.61 

2.19 

P . P 

1 .4 .'9 

1? 

1 7 . 4 

21 . 4 

208 

9 . 70 

2.32 

2.36 

.961 

?y . h 

21.2 

338 

15.95 

2.24 

P . •? 7 


13 

17.4 

X > , /■ 

203 

9 . .7 5 

2 . 2 It 

2 ■ 

1 .04 3 

29 . H 

21 . 4 

328 

15.33 

P . 1 _s 

2.18 

i £ ? 

14 

17.4 

21 .4 

202 

9.43 

? * pfi 

2.29 

1 . 126 

29.8 

21.3 

318 

14.94 

2.10 

2.12 

1 .898 

15 

17.4 

21 .6 

3 9i> 

9.00 

2.15 

2 . 19 

1.209 

29.8 

Pl.fi 

* t l ^ 

14.45 

2.03 

2 . 06 

2.05 

16 

17.4 

21 . 7 

19 B 

9 .03 

2.16 

2.19 

1 .292 

29.8 

21 .6 

308 

14.26 

1.99 

2 . 05 

2.18 

17 

17.4 

21 . 5 

187 

8.70 

2.08 

2.11 

1.373 

29.8 

21.fi 

307 

14.21 

1 . 99 

2 . 02 

2.32 

18 

17.4 


1 6:5 

P * f; 0 

2.06 

? . 0 9 

1 .455 

29.6 

PI .6 
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13.61 
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1 .94 

2.45 

19 

17.4 

2 1 . 7 

161 

B , It 1 i 

2.00 

2 .06 

I .538 

29.6 

21 .5 

289 

13.44 

1.69 
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2.59 

20 

17.4 

21.7 

1 63 

6.4? 

r* .01 

? .04 

1 -618 

29.6 

21 . 6 
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13.52 

1.91 

1.94 

2.71 

21 

17.4 

21 . 9 
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1.97 

2 . 00 
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29.5 

21.7 
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13 .04 

1 .84 
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2.84 

22 

17.4 

22.2 

171 
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1 .88 

1 .793 
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21.9 
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1.66 

1 . 71 

3.02 

23 

17.4 

22.2 
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1.68 

1 . 867 

29 . 7 
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12 .32 
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1.76 

3.14 
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TABLE II. - Concluded. EXPERIMENTAL DATA SUMMARY 
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Figure 1 

F n ( z ) 


.4 .6 

VELOCITY RATIO, z = u/u*, 

. - Hypergeometric functions for step-function analysis. 
= F(n-l+a, -n; a; z) ; a = (m+3)/(m+2); m = 5.6. 



















Figure 3, - Comparison of step temperature distribution data with analyses 




























CORRECTED MEAN Nl 




wall temperature (ref. 2). 















